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' In a recent significant advance, using Laguerre series, the valuation of Asian 

fN| , options has been reduced in [D] to computing the negative moments of Yor's 

q ■ accumulation processes for which functional recursion rules are given. Stress- 

or \ ing the role of Theta functions, this paper now solves these recursion rules 

Q . and expresses these negative moments as linear combinations of certain Theta 

integrals. Using the Jacobi transformation formula, very rapidly and very sta- 
bly convergent series for them are derived. In this way a computable series 
for Black-Scholes price of the Asian option results which is numerically illus- 
trated. Moreover, the Laguerre series approach of [D] is made rigorous, and 
vfH 1 , extensions and modifications are discussed. The key for this is the analysis 

J^J ■ of the integrability and growth properties of the Asia density in [Y], basic 

\ problems which seem to be addressed here for the first time. 

•5 

1. Introduction: Asian options are path-dependent options on the arithmetic average 
of the price of their underlying security. While they are widely traded financial securities 
their valuation still poses intriguing problems, even in the Black-Scholes setting. The 
Laguerre series approach of [D] so was a significant advance and reduced valuing Asian 
| options to computing the negative moments of the averaging process. For the latter it 
gave functional recurrence rules. Their structure and significance have been thoroughly 
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analyzed by Yor and his coworkers from a probabilistic point of view. 



Taking an analytic point of view, this paper explains how to express these negative mo- 
ments using Theta functions, and in this way derives computable Laguerre series for the 
value of the Asian option different from those of [SE]. The idea is that using the recur- 
rence relations of [D] the negative higher moments of the averaging processes should be 
linear combinations of certain Theta integrals. This is true in the basic case of §10, and 
depending on the relative constellation of risk neutral drift and squared volaltility only 
finitely many correction terms have to be added for the general case of §12. Using the 
Jacobi transformation formula for Theta series, §11 derives a series for these Theta inte- 
grals which is optimized with respect to speed and stability of convergence. We so obtain 
a Laguerre series for the value of the Asian option whose coefficients are series given by 
integration against Theta series. 

This occurence of Theta functions in valuing the Asian option is rather surprising. Indeed, 
higher moments of the averaging process have been thoroughly studied in [YE, Part A]. 
In their simplest form they were found to be expressed by Bougerol's identity as the ex- 
pectation of the respective powers of the hyperbolic sine evaluated on a Brownian motion. 
While a structural explanation of this link of Asian options with modular is still missing, 
it highlights a characteristic difficulty of the Laguerre series approach. Its workability 
crucially depends on the specifics of the option to be valued. 

For the Asian option, we think any valid expansion into orthogonal functions will even- 
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tually be based on Yor's Asia density of [Y, §6]. However, it is measurability which is 
addressed in [Y, §6] while for such purposes integrability is required. Thus we establish 
in §6 integrability of a certain class of functions of exponential type with respect to the 
Asia density, and put this result to work in deriving two Laguerre series for the value of 
the Asian option. 

First, we illustrate a natural principle for getting series expansions of option prices in terms 
of higher moments of the option's control variable. The idea is to Laguerre expand in the 
risk neutral expectation that gives this value the taking-the-non-negative-part function, 
and then try to put the expectation through the series so obtained. In §7 we combine our 
integrability results with characteristic mean convergence convergence results to justify 
this last operation, and so make this series available to valuing Asian options. 

Second, there are the Laguerre series in the spirit of [D]. Their idea is to construct the 
price of the Asian option as probability density in the strike price using the notion of ladder 
height densities, Laguerre expand this density, and try to re-interpret the Laguerre coef- 
ficients in terms of the reciprocal of the averaging processes. As explained in §8, Laguerre 
expandability of the Asian ladder height densities crucially depends on the integrability 
results of §6, and we have slightly modified the original series in negative moments of the 
averaging process of [D] in addition. 

Moreover, we give a second way of re-interpreting the Laguerre coefficients using expecta- 
tions of gamma function type integrands evaluated at reciprocals of the average process. 
These expectations also come up naturally if convergence of the Laguerre series in nega- 
tive moments is studied. Indeed, empirical evidence in §13 suggests that their convergence 
behaviour is controlled by a certain convergence parameter. Characterizing its optimal 
values will be in terms of these gamma type expectations. As indicated in §14, the latter 
now depend on Bessel functions, so have a different character than the negative moments, 
and we hope to return to their study elsewhere. 

2. Preliminaries on Laguerre series: This section collects pertinent properties of 
Laguerre polynomials from [L, §4] and [S] fixing any real number a > — 1. For any 
non-negative integer n, the a-Laguerre polynomial L%(z) is then given by 



ra/ 1 = V^ n ( ^k ^jn + a + l) z*_ 

n\ z ) 2^ k = { ] m+a+l) kUn- 



T(k + a + l) k\(n-k)\ 

for any complex number z. The first few a-Laguerre polynomials so are Lq(z) = 1, 
Lf(z) = 1+a-z, and L%(z) = (1/2) • ((l+a)(2+a) - 2(2+a)z + z 2 ), and for any positive 
intger n they satisfy the recurrence relation 



K +1 {z) = — L%(z) - ——K^z) 

n + 1 n+1 



The a-Laguerre polynomials are orthogonal on the positive real line with respect to the 
weight w a (x) = x a exp(— x) in that w a (x)(L c ^(x)) 2 L c ^ l (x) dx = T(n + a + l)/nl and 
Jo°° w a(x)L"(x)L ( ^ l (x) dx = for any non-negative integers n ^ m. To discuss when 
functions F on the positive real line have a Fourier-type expansion with respect to the a- 
Laguerre polynomials, consider the case where F is obtained by integrating up a function. 
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So suppose there is a function / on the positive real line which is integrable on any finite 
subinterval of the positive real line such that F(x) = F(0) + f f(y)dy, for any x > 0. 
Moreover assume that the functions ^w a F and y/w a +i f are square integrable on (0, oo). 
Then F is represented by the a-Laguerre series 



n=0 



which converges absolutely and uniformly for x in any subinterval [c, oo) of the positive 
real line. Any n-th generalized Fourier coefficient c n of this series is given by 



Alternatively, if F is continuous and y/w a F is square integrable on (0, oo), then F is also 
represented by the above Laguerre series, see [L, §4.23]. Convergence at zero of this series, 
however, needs a separate study, see for instance [S, p.366f]. 

3. Preliminaries on Theta functions: Following [M], the concept to be discussed is 
the classical Riemann Theta function •& given for any complex number z and any positive 
real number t by: 

This series converges absolutely, and uniformly on compact sets. Thus $ can be seen 
as a holomorphic function on the product of the complex plane with the upper complex 
half-plane. In modular forms this is usually done such that the above Theta function 
would be considered as evaluated not at t but at the point int of the upper half-plane. 
Theta functions are a basic class of holomorphic modular forms. They have been studied 
since the eightteenth century and are interrelated with a number of areas of central im- 
portance for mathematics, as number theory, algebraic geometry, classical analysis, and 
partial differential equations. 

Two such interrelations are to be described. For the first recall that t~ x l 2 exp (— nx 2 /t) 
is the fundamental solution of the Heat equation on the line with initial data at t = a 
delta function at x = 0. Thus ■& at real arguments can be seen as the superposition of 
infinitely many such solutions with initial data being delta functions at the half-integers 
x = n + 1/2. 

From the point of view of modular forms, Theta functions can be characterized by a cer- 
tain periodicity behaviour with respect to each of their two variables. Here, the behaviour 
with respect to the second variable is deeper and more subtle. It is expressed by a func- 
tional equation with respect to the second variable, the Jacobi transformation formula. 
Restricting to real arguments, this formula is the following remarkable identity: 

" X (-ire-(™) 2 -*cos(27mz) = -^y e -(*-W4 
=i y/irt ^nez 

for any real numbers x and t > 0. The left hand side of this identity, the Jacobi transform 
of is rapidly converging for t large, whereas its right hand side is rapidly converging for 
t small. The Jacobi Transformation Formula is proved by Fourier analysis, exhibiting $ as 
the Fourier expansion with respect to the second variable of its Jacobi transform. 
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4. Basic notions: First we discuss the notions basic for the analysis of Asian options. 
We work in the Black-Scholes framework using the risk-neutral approach to the valuation 
of contingent claims. In this set-up there are two securities. There is a riskless security, a 
bond, whose price grows at the continuously compounding positive interest rate r. There 
is also a risky security, whose price process S is modelled as follows. Consider a complete 
probability space equipped with the standard filtration of a standard Brownian motion on 
the time set [0, oo). On this filtered space, we have the risk neutral measure Q, which is a 
probability measure equivalent to the given one. And then we have a standard Q-Brownian 
motion B such that S is the strong solution of the following stochastic differential equation: 

dS t = w • S t • dt + a • St • dB t , t e [0, oo). 

The positive constant a is the volatility of S. The specific form of the otherwise arbi- 
trary constant w depends on the nature of the security modelled (eg. stock, currency, 
commodity etc.). For example it is the interest rate if S is a non-dividend-paying stock. 

Fix any time to and consider the accumulation process J given for any time t by: 



J(t) = S u du , 

Jtn 



The European-style arithmetic-average Asian option written at time to, with maturity T, 
and fixed-strike price K is then the contingent claim on the closed time interval from to to 
T paying ( J(T)/(T-t ) - K)+ := max{0, J(T) / (T-t ) - K} at time T. Recall that points 
in time are taken to be non-negative real numbers. The price Ct of the Asian option at 
any time t between to and T is given as the following risk neutral expectation 



C t = e- r <> T - t )E < t 



T-t 



which is conditional on the information Tt available at time t. However, following [GY, 
§3.2], do not focus on this price. As described there in great detail, we instead normalize 
the valuation problem, consider the factorization: 

-r(T-t) 4 a 

C t = - =±-CM(h,q), 

T-t a 2 V lH> 

and so reduce to computing 

C^(h,q) = E^[(A^-q)], 

the normalized time-t price of the Asian option. To explain the notation, is Yor's 
twohundred percent volatility accumulation process 

A W) = f h e l(B w +»w) dWj 

Jo 

and the normalized parameters are as follows: 

u=^-l, h = ^(T-t), q = kh+q\ 
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where 



k 



K 
~St 



q* = Q*(t) = ^(k- (t-t ) - J* S u dv)j 



To interpret these quantities, v is the normalized adjusted interest rate, h is the normalized 
time to maturity, which is non-negative, and q is the normalized strike price. 



5. A reduction of the valuation problem: Computing the normalized time-t price 
of the Asian option reduces to the case where the normalized strike price q is positive. 
Indeed, if q is non-positive, Asian options loose their option feature, and their normalized 
time-t price is given by 



C (u) :=E Q[( A M =E Q 



-Q- 



On applying Fubini's theorem, this last expectation is computed as follows: 



E Q 



A 



(")' 



e 



2h(v+l) 



1 



2{v+\) 



for v any real number. The right hand side is analytic in v with its value at v = — 1 being 
equal to h. For computing C(") if q is positive we discuss two Laguerre series. They in 
turn crucially depend on an extension of a basic result of Yor's about the Asia density to 
be discussed next. 



6. Yor's Asian option density revisited: This section discusses growth and integra- 
bility properties of the density of the accumulation processes A( v \ This is based on Yor's 
integral representations of [Y, (6.e), p. 528] recalled below and extends his measurability 
results. Recall the Girsanov identity for the normalized time-t price 

C<"> = EQ [(Atf - q) + ] = e~^E[(A^ - q)+e uWh ] 

which follows on making Wt = ut + Bt a Brownian motion and dropping reference to the 
resulting measure. It reduces computing CM to a corresponding problem in the zero 
drift accumulation process at the cost of introducing a second stochastic factor in the 
integrand. Define the density a Vi h on the positive real line by 

roo xy 2 

a Vih (x) = e 2 / g v {y)e 2 tp xy (h) dy , 
Jo 

for any x > 0. Here g v is the map on the positive real line that sends any y > to its u-th. 
power, and ip^(h) is for any positive real number ^ given by: 

ip%(h) = J e _2 7T ^ cosh ^ w \inh(w) sin (J^-^j dw . 

To address integrability let / be any continuous function on the positive real line, and for 
any real numbers a, b define the weighted maps f a ^ by f a ,b( x ) = ex P (x~ 1 a)x b f(x), for 
any positive real number x. Then we have the 
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Proposition: If a < 1/2 and f behaves like a power map near the origin and towards 
infinity, the map sending any positive real number x to f(x~ x ) exp (ax)x a v ,h{x) extends 
to an integrable map on the non-negative real line, and we have the identity 

/■OO -1 

E[f a , b (A^)g v (e Wh )] = c h / f a , b (-)a v , h (x) dx 

Jo Kx/ 

of finite integrals, where Ch = (27r 2 /i) -1 / 2 exp ((2/i) _1 7r 2 ) . Moreover, the map sending any 
positive real number x to exp (ax)x~ b a Uj h(x) extends to an integrable and bounded map on 
the non-negative real line. 

Remark: Examples of functions / which meet the conditions of the Proposition are fur- 
nished by any rational function with no poles on the positive real line, and by any function 
sending x > to (x — q) + or (q — x) + . 

The proof of the Proposition is based on Yor's triple integral [Y, (6.e), p. 528] and his 
observation [Y, (6.g), p. 529]. Recall the latter implies that for any non-negative integer k 
the map £ i— > ip$(h) is of order big Oh of £ fc with going £ to zero, while the former asserts 

/•OO /"OO 1 X y 2 

E[h(A^)g(e Wh )} = c h / / M-Je^ ^(y)e~2" ^ xy (h) dx dy 

Jo Jo Kx/ 

for any Borel measurable functions g, h on the non-negative real line into itself. In the Yor 
triple integral for the Proposition you would like to interchange the order of integration. 
With the functions / and £ i— > ^>^(h) possibly negative, Tonelli's theorem does not apply 
for this. However it does apply on taking absolute values of the integrands. And if the 
resulting integral can be proved finite this will give integrability and justify interchang- 
ing the order of integration in the original Yor's triple integral using Fubini's theorem. 
Abbreviating 7 = v — b and 5 = 1/2 — a we are so reduced to prove finiteness of 

POO / \ \ 

J = J \f\( K -je~ Sx D 0jOO (x)dx, 

setting 

b ?2 

Da,b{x)= [ e -^-C-\Mh)\dC. 

J A x 

Decompose J into four parts by breaking the two integrations at the point 1. The idea is to 
first majorize the respective inner integrals Da,b using Yor's above observation, and then 
use the exponential decay of the integrand towards infinity to bound the whole summand. 
We discuss this for the case where x and £ both range from 1 to infinity. Fixing x > 1 and 
changing variables n = (2x)~ 1 ^ 2 r] : we have 

x ' 

To majorize this last integral notice that (2x) 1 / 2 r] is greater than or equal to 1 by con- 
struction. On [l,oo) the map y 1— > y 7 exp(— y) is bounded by a constant A 1 depending 
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only on 7. Use this in the defining integral for tp to obtain 

f°° _w 2 _ f°° _u£_ 

y 1 \ip y (h)\ < I e 2h y 1 e y sinh(w) dw < A 7 / e 2h sinh(w) dw . 
Jo Jo 

Thus there is a constant B 1)h depending only on 7 and h such that D lj00 (x) sits below 
B-y t h ■ x~^ +1 ^ 2 ^ for any x in [1, 00). For the respective summand of J we so obtain, 



/oo -1 poo -1 "2 (5 

1/1 (-)e- Sx D hoo (x) dx < B,, h ^ l/l (-) ■ A_ . . 

Now apply the above reasoning to the middle factor of the last integral's integrand and 
bound it by (5/2) 7+1 / 2 A_( 7+1 / 2 ). On substitution, 

/OO -I 
|/|(i)e- fa I? 1>00 (x)dx 

/£n 7+1/2 /•<» /U _«5 T 

<5 7 ,.-(-) A_ (7+1/2) ^ |/|(-)-e ^cfa, 

and this last integral is finite since 5 is positive by construction and / behaves like a power 
map near the origin by hypothesis. The remaining three summands of J can be majorized 
in a similar way, and this completes the proof. 

7. First Laguerre series for the Asian option: As a first application of our study 
of the growth behaviour of Yor's Asia density we illustrate a general priciple for getting 
series expansions of option prices in terms of higher moments of the option's control vari- 
able. The idea is to Laguerre expand in the risk neutral expectation that gives this value 
the taking-the-non-negative-part function, and then try to put the expectation through 
the series so obtained. The krux of this Ansatz is to justify the last operation, and if such 
a justification can be found it will be specific to the option considered. For the Asian op- 
tion we have been able to give such a justification using the properties of the Asia density 
established in §6 Proposition in an essential way. The precise result is the following series 
in terms of the moments of the reciprocal of Yor's accumulation process. 

Lemma: Let c be any positive real number with < qc < 1/2. Setting X = A h , the 
normalized time-t price C(^) of the Asian option has the absolutely convergent series 



n=0 

whose any n-th coefficient c n is given by 

\k+l 



^ = E rj^I) Q ((*+«> •7(* +Q+ i,c) + ^-V=), 

with 7(s, a) = exp (— x)x s ~ 1 dx for any a, s > the incomplete gamma function, and 
E Q 



X - L Ax)\ } T(k+a + l) k\{n-k)\ E 



X k-i\ 
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Remark: In [D, Theorem 5.1, p.417f] a Laguerre series for the the density of the recipro- 
cal of Yor's accumulation process is given. Use it to compute the expectation defining 
by a formal term by term integration. The series so obtained is then seen to be identical 
with that of the Lemma. 

The proof of the Lemma illustrates the basic principles for working with orthogonal poly- 
nomials. Take any function 4> on the positive real line which has an ct-Laguerre series, so 
let w x l 2 § be square integrable where w is the weight w(x) = exp(— x)x a . From [AAR, 
Theorem 6.5.3, p. 308] we have mean convergence with respect to w of this series to <p. 
Writing Ck for any k-th Laguerre coefficient and R n = 4> — (cqLq + ■ ■ ■ + c n L^) for any 
n-th order remainder term of the Laguerre series, this is equivalent to 



(•OO 

lim / w(x) ■ [R n (x)] 2 dx = 



The idea for the proving the Lemma is to transcribe the ordinary convergence of the series 
to be proved into mean convergence of the Laguerre series for 4> c (x) = (c— x) + and apply this 
mean convergence result. Thus let R n now denote any n-th order remainder term of the 
Laguerre series for (j) c and set p n (x) = xR n (x~ 1 qc). Observing (x — q) + = c~ 1 xc/) c (x~ 1 qc). 



<k=0 



E Q 



qc s 
X, 



c 



The proof of the Lemma so reduces to show that this last expectation goes to zero with n 
to infinity. Using the Asian density a v ^ of §6 Proposition, we have 



E Q 



Pn(X) < C J a^ h (x)p n (^j 



dx 



with Co = Ch exp (— u 2 h/2). Applying Cauchy-Schwarz, 

"OO 



E Q 



Pn(X) 



<4 







1 



ar 



a„ jh (x)dx- / \R n (qcx)\ a Ujh (x) dx 



o 



/•oo / C \ 

<'^[X-'}- \RJS)\ 2 a„, h (1 )«, 
qc j \qc/ 



on changing variables ^ = qcx for the last equality. Since all higher moments of any Yor's 
accumulation process are finite using [Y, §4], we are further reduced to show that the 
second integral factor goes to zero with n to infinity. This is where the growth proper- 
ties of the Asian option density enter in an essential way. Indeed, using the hypothesis 
< qc < 1/2 choose any qc < (3 < 1/2. Then §6 Proposition applies to show that 
x i— > exp {(3x)x~ a a v ^h{%) is bounded on the positive real line, whence 

is bounded on the non-negative real line as well. As a consequence there is a positive 
constant D such that 

[°° \R n (0\ 2 a„J^)dZ<D e-H a \R n (0\ 2 d^ = D f°° «;(£)|i? n (£)| 2 d£ . 
Jo X Q C/ Jo Jo 
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Now the mean convergence result recalled at the beginning applies to give that this last 
integral goes to zero with n to infinity. Absolute convergence of the series follows using Rie- 
mann's criterion for conditional convergence: with the above series also any rearragement 
of it converges to the same function. This completes the proof of the Lemma. 

8. Ladder height density Laguerre series for the Asian option: A second La- 
guerre expansion of the value of the Asian option has been proposed in [D] . The idea is to 
construct this value as a probability density in the strike price, Laguerre expand, and try 
to re-interpret the Laguerre coefficients in terms of the reciprocal of Yor's accumulation 
process. The key notion is that of the ladder height density g z associated to a probability 
density function gz of any positive Q-integrable random variable Z. It is defined for any 
c > by 

1 



9z ^ = EQ[Z]J w -9z{w)dw. 

We have slightly modified the argument of [D, §7] how this applies to Asian options. Using 
§5 assume the modified strike price q positive. For any positive real number c, define the 
Q-probability density function g of a positive and Q-integrable random variable by 

g(y) = r , ■ where Y = 



E[Y~i] y ~ A W)' 

Then we have the representation of the normalized time-t price of the Asian option 
in terms of the ladder height density g associated to g 

C (u) =E Q^_ q + t^ E Q . = (c) 

setting X = Postponing the computation of .EfX -1 ] to §12 Lemma and with E®[X] 

computed in §5, we are so reduced to computing g(c). The idea is to try to use a Laguerre 
series for this. Thus consider for any positive real number c and any real numbers (3 and 
5 the formal Laguerre series expansion 

c P e~ 5c ■ f(c) = J2^ =0 c n L*(c) 
whose any n-th Laguerre series coefficient is given by 

(-l) k (r? 



En 
k=0 



T(k+a+l) \k 
where I/3,s,k abbreviates the integral 



lf3,5,k 



T / a+/3+k -(1+5)2/ =/ \ 7 

Ip,s,k = V e g(y) dy . 

Jo 



The basic question is under which conditions on (3 and 5 this series converges, and this is 
again crucially based on the analysis of Yor's density in §6. The second question is how to 
compute the coefficient integrals Ip : s,k for these admissible values of j3 and 5. Here is the 
following extension of [D, Theorem 7.1, p. 422]. 
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Proposition: Assume q > 0, abbreviate X = A h , and let f3 and 5 be such that 
a + 2(3 > — 1 and 26 > — (1 + (2qc)~ 1 ). Then the above a-Laguerre series converges. 
If moreover qc < 1/2, we can choose 5 = — 1 as an admissible such value, and 

1 

Xa+P+k+l ' 

for any non-negative integer k. If 5 ^ — 1 and m = a+(3+k is any non-negative integer, 



_ 2 (qc) a+ P +k 
~ EQIX- 1 } ' (a+P+k + l)(a+P + k+2) 



E Q 



Iqc 



mi 



m—£ 

-Y 



mi 



1 1=0 

(qc) m - i - p 



E[X~i] 



(1+5Y+ 1 {(l+5)™- £ 
1 



-p)l(l+5)p+ 1 



E 



Xm—t—p—l 



qc(l+S). 

~x — 



To prove the Proposition, given [D, Theorem 5.1, p. 418] we first have to check square 
integrability, i.e., have to determine which maps x i— > x a+2/3 exp ( — (1 + 25)x)(g(x)) 2 are 
integrable on (0, oo). To get a hold on g, recall if Xj is any random variable with Q 
density function g, we have E®[X-g]g(x) = Q( exp (aXg) > exp (ax)) for any a > 0. Hence 
Markov's inequality gives E®[Xg]g(x) < exp (— ax)E®[ex.p (aXj)]. On substitution we 
are so reduced to characterize finiteness of 

E Q [exp (aXg)] f°° ^ a +2/3 e -(i+2(s+a))x dx 



The integral factor gives the conditions a + 2/3 > —1 and l + 2(a + 5) > 0. Since 
E Q [Xj] = E Q [X 2 ]/2 using [D, Theorem 2.1, p. 411] and all moments of X are finite, 
we have to determine those a > for which E^[exp (aXj)] is finite. Since on unravelling 
definitions 



£0[exp(a%)1 ^(_i_ £ « 



Xe 



aqc 



- 1 



this again is an application of §6 Proposition. The expectation will be finite for any posi- 
tive real number a such that aqc < 1/2 yielding 25 > — (l + (2qc)~ 1 ). This is satisfied for 
5 = —1 if qc < 1/2, and in this case the coefficient integrals Ip,s,k have been computed in 
[D, Theorem 2.1, p. 411]. An analogous computation gives the coefficient integrals in the 
other case which is not contained in [D]. This completes the proof of the Proposition. 



9. Review of two results of Dufresne: As one of his insights into the structure of 
the accumulation process 

AM= f e^+^dw, 
Jo 

Yor describes in [Y, (4.k), p. 522] how its higher moments are determined by certain Gauss 
hypergeometric functions 2-Fi- These results have been extended in [D, §4] to the higher 
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moments of the reciprocals of these accumulation processes in the following way. First we 
have the integral representation 



E Q 



1 



4"> 



2 - / ye 27T v \ - ,y > dy, 



for any real number z/. At the base of this result is Yor's description of the Laplace trans- 
form with respect to time of the first moment of A( v ) using the confluent hypergeometric 
function $. This description remains valid mutatis mutandis for the expectation of r-th 
powers of A( u ) where r > — 1. The identity resulting on inversion is in terms of the Gauss 
hypergeometric function and holds for r > —3/2 using analytic continuation. 

The higher moments of the reciprocal of A( u ~) are determined recursively from its first 
moment. Considering any k-th such moment 

m k (h) = E^[(A^)- k ] 
as a function in the time variable h > 0, we have 

m k = 2(k - (i/+l))m fc _i - -j^j m 'k-i » 

for any k > 2. This is proved using the Ito Lemma on applying time reversal. Different 
angles on this are discussed in [Yl], [Y2], and [Y3] in particular. 

10. Higher moments as Theta integrals in the basic case: It is Theta functions 
as discussed in §3 which provide the key for understanding the higher moments of the recip- 
rocal of Yor's accumulation process. With higher moments of Yor's accumulation processes 
themselves being determined by hyperbolic sines [Y,4.g, p. 519], this appearance of theta 
functions seems unexpected at least. This section considers the basic case \v — 1| < 1 in 
which the higher moments of the reciprocal of any time-/i value of A( u ) 

mk (h) = EQ[(A^y k } 

are completely determined by theta functions. Indeed, computing them reduces to com- 
puting the Theta integrals <d k given by 

1 f°° (v \ w k ~ x 
1 ' 2V2 Jo V2 J (wh + l/2) k + 1 / 2 

for any h > and for any positive integer k. The precise result is the following 
Lemma: Suppose \v — 1| < 1. For any positive integer n, we have 

e^~m n (h) = a n ,i@i(h) H h a n , n Q n {h) , 
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for any h > 0, where the coefficients a Uj k are recursively determined by a\^\ = 1 and the 
recurrence relations 



dn+1,1 — 
0>n+l,n+l = 




for any k between 2 and n. 

Remark: Similar but different charactarizations of these higher moments have been ob- 
tained by Yor and coworkers if v = 0, see in particular [Yl], [Y2], and [Y3, §4]. 

The key insight for the proof is that in Dufresne's integral for the expectation recalled in 
§9 the hyperbolic quotient factor of the integrand is the Laplace transform at y 2 of a theta 
function if \v — 1| < 1. Indeed, for any positive real number y we then have 



cosh((z/ — l)y) 
ysinh(y) 




Applying Fubini's theorem the expression of the Lemma for m\{h) follows from this. An 
induction using the recursion relation for mt of §9 then completes the proof. 



11. Computing the Theta integrals: This section discusses two series for computing 
the Theta integrals © n . From §10 these are recalled to be given by 

1 f°° (v \ w n ~ x 

©n(M = 7= / $[—w)- — —T^rdw, 

1 ' 2V2 Jo ^2 J (wh + l/2)»+V2 ' 
for any h > and for any positive integer n. 

First series: First we have the following absolutely convergent series in terms of the 
confluent hypergeometric function of the second kind \& 

^ /,n r(n) T / 1 ^\ r(n) v^oo / 1 (7rm) 2 \ 

Q n (h) = — In, -; + > cos(7rmz/) • * In, -; - — r— . 

v ' 2h n V ' 2 ' / /i n ^m=i v ; V ' 2 2h J 

It is obtained by term by term integration of the series which results from applying the 
Jacobi transformation formula of §3 to the Theta function factor in Q n (h). However, it 
is problematic in that confluent hypergeometric functions are most complicated functions 
and their implementations are often flawed. Similar series for v = are in [Y3, §4]. 

Second series: Still there is a second series in terms of less problematic functions. It is 
the absolutely convergent series 

00 

®n(h) = ^ C B,m(h) + ^ d B,m(h) 
m=0 m£Z 
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whose coefficients are for any positive real number B given as follows: 

n-l / 



1 



-1 



fc=0 n 
n-l 



c B , m (/i) = 7m £(-2) 



fc=0 



2 

71 — 1 



n — 1 



1 



1 V k J (2/i J B+l)"- fc - 1 /2' 



fe / (7rm) 2fc 



^(1) ^(2) 



for any positive intgeger m, abbreviating 

cos(nmu) ( — l) n_1 



,2n-l 



V2 



2 n-l /j2n-l/2 



c. 



(1) 

m,k 



-k-1 



-(7rm) 2 B 



(2) 
m,k 



•1) 



=0 
n—k 



(4). 



^in-fc ( (7rm) 2 ( S + 275; 



1 ' 



-{irm) 2 B 



w(nm^B+ (2K)- 1 



If m is any integer different from —v/2, we have 



(2) 
m 



abbreviating 



(i) 

m,k 



m + - 



2(n-fc-l) 



-i(-+f) ; 



(J) 



(2)n (| + 2/i) 



fc+1/2 



l) n Z/ 2n-l 

mH — 
2 



v 

m + - 



v / B~ 1 + 2h , 



while for m = — 1//2 any integer, we have 

2 n-l 



^n-l/2 



0F(n-l) (2^+l)-V2 



Here the function W is given by W(z) = exp (z 2 )Erfc (2), for any complex number z, and 
(A)fc is the Pochhammer symbol of any complex number A recursively given by (A)o = 1 
and (A)fc_|_i = (A — k) ■ (A)& for any non-negative integer k. The idea for this series is to 
optimize the convergence behaviour using the Jacobi transformation formula of §3. It is 
thus obtained by breaking the integral defining Q n (h) at B and integrating over [B, 00) 
now the Jacobi transform of the Theta function in Q n (h). 
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12. Higher moments in the general case: Computing the moments of the recipro- 
cal of any time-/i value of A( u ) 



m k (h) = E[(A^y k ] 



for arbitrary indices v is not solely in terms of the Theta integrals Q k (h) of §11. Additional 
non-Theta correction terms are required. The precise result is the following 

Lemma: Let n* be the smallest integer m such that 2m+l— is non-negative. For 

any positive integer N , we then have 



e 2 m N (h) =J~] a Nyk e k (h) 

< ' k=l 

En* — 1 ^-^7V— 1 r -\ 
n=0 ^k=0 { bN M h )- C nM h )- a N,k+l-Dn,k+l(h)j, 



for any h > 0. 

The functions C n ^ k and D Ujk : The functions C n ^ k and D njk are for any positive real 
number h given as follows. First, 



C n , k {h) = h*(^ 



2k+l 2fc+l 



Here Wg is given in terms of the complementary incomplete gamma function T(a, x) by 
Wg(x) = exp(x 2 ) • r((£+l)/2, x 2 ) for any real number x, and (3 n = 2n + l— — Then, 



D n>k (h) = 2 



-1 (~27n)' 
(\)k 



^(v^T)-fg(i) 



2 • (h ln y+i 



where W is the function Wo, abbreviating 7 n = 4 _1 (|z/ — 1| — 1 — 2n) 2 , and with (A)fc 
the Pochhammer symbol of any complex number A recursively given by (A)o = 1 and 
(A)fc+i = (A — k) • (X) k for any non-negative integer k. 

The coefficients a n ^ k and b nyk : The coefficients a n ^ k are those of §10 Lemma. The 
coefficient functions b U}k (h) have the form 



l /M-A ,-(§+2(n-l)) 

O n ,n-l{ n ) — °n,n-l,n-l" 

and for any non-negative integer k less than or equal to n — 2 



n-l 



-(l+fe+£) 
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They are recursively defined by 



h,o(h) = 



2 -s 
h 2 . 



b 2 , (h) = ^(2(l-v) + ^-)h- 1 + 



3 h~* 



and t>2,i(h) 



2 -i 
h \ 



and the following coefficient recurrence rules for n > 2 



bn,n—l,n—l i 

n 



( 



i,n-i,e 



f[b njn -2,£-l 
,2 



V \ 1 

2(?1 — Z/) + 2^ j^n,n-l,n-l — J^>n,n-1,n-\ 

> ^(|+ 2 ( n_1 )) & n,n-l,n-l ~ ^K,n-2,n-l 



For any positive integer less than or equal to n — 2 we have 

b n +i,k,e 

(2(n-u) + ^jb n , k ,i+ ±(§+fc+*-l)& n , fc) /_i - ^n,n-i, 



£ = 

£ = l,...,n-2 
£ = n-1 
£ = n. 



= < 



e-i 



and 



b n +i,o,e — < 



£ = 

1 = 1,. 
£ = n, 



' (2(n-v) + ^y n , 0fi £ = 

(2(n-u) + ^jb nAe + 1 (%+£-iy n ,o,£-i e = l,...,n-l 



. , n — 1 



13. An example: This section considers the example of valuing an Asian option on a 
non-dividend-paying stock with an interest rate za = r = 9% p. a., a volatility a = 30% 
p. a., and with time to maturity 1 year from today time 0. Suppose today's price of the 
stock S is equal to 100, and the option has been issued today at par, i.e., K = S . Then 
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v = 1.0, h = q = 0.0225, and we have for the Theta integrals Q n (h) of §10: 



n 


n c ,30 


nd,30 


e n (/i) 


1 


4 


4 


44.2790749547 


3 


4 


4 


46822.1151330265 


5 


4 


4 


70467169.1745463509 


7 


4 


4 


116806496442.2335048967 


9 


4 


4 


202679488380050.1455046663 


11 


4 


4 


361222145817967427.8404500431 


13 


4 


4 


655202294765472642150.1383863999 


15 


4 


4 


1203316829766524026893611.5309617920 


17 


4 


4 


2230542763985591766299300393.1676737371 


19 


4 


4 


4164445349148560860975785222305.1797500863 



Table 1 Growth of Theta integrals Q n (h). 



We have computed these values using the second series of §11 with B = 0.3, and denote 
by n Cj so and nd,3o those indices m for which the first 30 decimal places of the partial sums 
c s,o + ' ' ' + CB jm and ds,o + - • - + dB,m respectively are correct. The calculations have been 
done on a HP Visualize' C200 using the GP/PARI CALCULATOR Version 2.0.11 (beta) 
of C. Batut, K. Belabas, D. Bernardi, H. Cohen and M. Olivier with 150D precision. Such 
a precision was desired since the coefficients a U} k of §10 Lemma grow like k factorial. Using 
this result we compute the following values of the n-th moments m n (h) of the reciprocal 
of Yor's accumulation process A( u ~) at time h: 



n 


m n (h) 


1 


43.7837269185 


3 


91741.4898743896 


5 


215689033.4795106594 


7 


566731384819.7874280500 


9 


1657943864684789.3944805424 


11 


5380692663910949427.6422561855 


13 


19305805471114617878830.4144899460 


15 


76329904185193047892084144.9033200670 


17 


331513545523453183819373799801.9730752657 


19 


1576936784103137595901225118618858.9370371494 



Table 2 Growth of the moments m n (h). 



Given this rapid growth of the negative moments the question is if and when convergence 
of the Laguerre series is rapid enough so that only the first few of their terms are needed 
to reproduce the price of the Asian option. Principally, convergence of Laguerre series 
is effected by the oscillation of the Laguerre polynomials in the range of positive real 
numbers smaller than 4n+ 2(ct+l). In both of our Laguerre series this is controlled by the 
parameter c in the sense that the smaller c the more of their terms are in their respective 
oscillatory range, and the better should be convergence. The Black-Scholes price of our 
Asian option is given as C BS = 8.83 in [RS, Table 1.3] from which its normalized time-t 
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price Ct") = 0.002173850758 is obtained on divison by 4061.916379. In the sequel we 
consider the ladder height density series of §8 Proposition with 6 = — 1 and (3 = and take 
a = 0, i.e., compute with the classical Laguerre polynomials. Using the obvious notation, 
we record in the following examples the values the series give for these two prices after 
summing the first n+1 terms, and give the respective amounts A n that have to be added 
to obtain their above target values. First we choose for c a "reasonable" value, to improve 
the series' sprint qualities so to speak. 



n 


r (y) 




r<Bs 


A BS 


1 


0.0041430388 


0.0019691881 


16.82867729 


7.9986772963 


3 


0.0029964651 


0.0008226144 


12.17139078 


3.3413907866 


5 


0.0030265291 


0.0008526784 


12.29350849 


3.4635084905 


7 


0.0027607453 


0.0005868946 


11.21391670 


2.3839167088 


9 


0.0024684544 


0.0002946036 


10.02665531 


1.1966553198 


11 


0.0023864317 


0.0002125809 


9.69348587 


0.8634858736 


13 


0.0023930119 


0.0002191611 


9.72021406 


0.8902140555 


15 


0.0023544102 


0.0001805595 


9.56341738 


0.7334173791 


17 


0.0022627357 


0.0000888849 


9.19104301 


0.3610430106 


19 


0.0021738508 


4.049 x 10" 15 


8.83000000 


1.645 x 10" 11 



Table 3 Ladder height series with c = 1.367054258545 



As to be expected for this magnitude of c, there is a steady convergence of the series. One 
should expect it to converge mildly oscillating. However, mind the following fallacy. It is 
not sufficient to have the the normalized time-t price CM correct up to the the first two 
or three decimal places only! Indeed, we have this accuracy for C^) from the outset while 
there are still big deviations from the Black-Scholes price C BS . These phenomenons are 
much more pointed if a larger value for c is chosen. Indeed notice how the following series 
actually start with almost the correct values, then wildly oscillate, and only begin to calm 
down after fifteen terms or so. 



n 


r M 


A M 




a bs 


1 


0.0020571929 


-0.0001166579 


8.35614534 


-0.473854662 


3 


0.0045737040 


0.0023998531 


18.57800280 


9.748002806 


5 


0.0044363733 


0.0022625225 


18.02017734 


9.190177342 


7 


-0.0055198391 


-0.0076936899 


-22.42112502 


-31.251125020 


9 


0.0028067490 


0.0006328989 


11.40077951 


2.570779513 


11 


0.0034914141 


0.0013175634 


14.18183220 


5.351832204 


13 


0.0017503210 


-0.0004235298 


7.10965752 


-1.720342482 


15 


0.0015816389 


-0.0005922118 


6.42448511 


-2.405514894 


17 


0.0017276113 


-0.0004462395 


7.01741241 


-1.812587586 


19 


0.0017603640 


-0.0004134868 


7.15045128 


-1.679548715 



Table 4 Ladder height series with c = 6 



In summary, we are now able to efficiently compute the negative moments of Yor's accu- 
mulation processes that enter into the Laguerre series for the Asian option. However, there 
has still work to be done in characterizing optimal choices of the convergence parameter c. 
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14. Epilogue: In this paper we have concentrated on the first of the two ladder height 
density Laguerre series of §8 Proposition. Working with the second of these requires to have 
for certain positive real numbers a, (3 at least computable expressions for the expectations 



faAt) = E Q 



X~ a e 2X 



for any time t > setting X = Af\ Adapting the argument of [Y, §3, p.515f], their 
Laplace transform with respect to time is computed in terms of modified Bessel functions 



In and MacDonald functions as follows 



2 Z" 00 

£(/«./0(*) = / C-'l^OKa-uivO dZ 

V Jo 

setting 7] = ((l+/3)/2) 1 / 2 and with z any complex number whose real part is positive and 
such that Re((2^) 1 / 2 ) > 1 — v. In analogy to the modified Weber-Schafheitlin integrals 
of [W, p. 410] these Laplace transforms can be expressed using the Gauss hypergeometric 
function 2 F\ • On Laplace inversion f a ^ is seen to be given by the contour integral 

' V a ~ u 2™ Jio S c R V2^W Jo 
Here log Cr is the following path of integration 

log R+in 





1 > 


logi? 


l s- 

1 < 



log R—iir 



Figure 1 The contour logC H 



where the positive real number R has to be chosen so big that the real parts of w 2 and 
w of any of its elements w are positive and such that rj is bigger than cosh(log-R). Under 
additional hypotheses, from [E, 7.7.3 (26), p. 50] the inner integral can be expressed using 
the Gauss hypergeometric function 2 F\. Alternative angles on these expectations, like 
their representation as Kantorovich-Lebedev transforms, are discussed in [AG] . 
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